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ABSTRACT 

The formation of quasar broad-line region (BLR) clouds via thermal instability in 
the presence of Alfven heating has been discussed by Gongalves et al. (1993a, 1996). 
In particular, these studies showed the relevance of Alfven heating in establishing 
the stability of BLR clouds in the intercloud medium. The present paper shows the 
results of time-dependent calculations (we use a time-dependent hydrodynamic code) 
following the evolution of BLR clouds, since their formation from the 10 7 K intercloud 
medium. We also calculate the UV and optical line emission associated with the clouds 
in order to compare with observations. Our results are compared with those of UV 
and optical monitoring of well-studied AGN, which suggest that the BLR is most 
probably composed of at least two different regions, each one giving rise to a kind of line 
variability, since low- and high-ionization lines present different patterns of variability. 
We discuss the alternative scenario in which lines of different ionization could be 
formed at the same place but heated/excited by distinct mechanisms, considering as 
a non-radiative mechanism the Alfven heating. 

Key words: Quasars: BLR - Hydrodynamics - Plasma: Alfven heating - Emission- 
line variability 



1 INTRODUCTION 

The study of the broad-line regions (BLR) of active galac- 
tic nuclei (AGN) is crucial in understanding the central en- 
gine of these objects. In addition, as the BLR reprocesses 
energy emitted by the continuum source at ionizing ultravi- 
olet energies that cannot be observed directly, the monitor- 
ing of emission lines provides unique pieces of information 
on the central continuum source. Photoionization by con- 
tinuum radiation from the central source is almost certainly 
the most important heating/ionizing mechanism in BLRs. 
One should note that non-thermal processes, for instance 
incoherent synchrotron radiation, were earlier invoked to ex- 
plain AGN spectra (for instance, Rees 1984, 1987; Contini 
& Viegas-Aldrovandi 1990). In any case, the pure photoion- 
ization model has a number of inconsistencies — the 'energy 
budget problem', the 'line ratio problem' and the 'line vari- 
ation problem' (Dumont, Collin-Souffrin & Nazarova 1998). 
In fact, Dumont et al. (1998) speculate on the need and na- 
ture of a non-radiative heating source of energy — dissipa- 
tion of sound or Alfven waves, weak shocks, magnetic dissi- 
pation, etc., as a possibility to partially solve these problems. 
The non-radiative process discussed here is the interaction 
between Alfven waves, which results in a dissipation pro- 
cess heating the gas. The present study aims at evaluating 



the effects of Alfvenic heating (AH) as non-radiative heat- 
ing, in addition to pure photoionization, on the broad-line 
spectrum. 

We studied two damping mechanisms for Alfven waves: 
the nonlinear and turbulent. These damping mechanisms 
have been investigated before, showing the relevance of 
Alfven waves and their dissipation processes in all the 
following astrophysical objects. Late-type and proto stars 
(Jatenco-Pereira & Opher 1989a, b; Vasconcelos, Jatenco- 
Pereira & Opher 2000); solar wind (Jatenco-Pereira & Opher 
1989c; Jatenco-Pereira, Opher & Yamamoto 1994); galactic 
and extragalactic jets (Opher & Pereira 1986; Gongalves, 
Jatenco-Pereira & Opher 1993b); early- type stars (dos San- 
tos, Jatenco-Pereira & Opher 1993; Gongalves, Jatenco- 
Pereira & Opher 1998); broad line regions of quasars 
(Gongalves, Jatenco-Pereira & Opher 1993a, 1996); and fi- 
nally, optical filaments in cooling flows (Friaga et al. 1997). 

The scenario for BLRs considered in this article is the 
two-phase BLR, in which the optical-UV lines are emitted 
by clouds embedded in a hot, external, intercloud medium 
(Krolik, McKee & Tarter 1981). However, thermal pressure 
alone is not enough to confine BLR clouds if the exter- 
nal medium is heated only by Compton scattering (Math- 
ews & Ferland 1987; Mathews & Doane 1990). Mathews 
& Doane (1990) take into consideration various physical 
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processes in the heating-cooling the BLR, including: i) re- 
combination heating and cooling; ii) radiative losses in the 
range between 10 4 K to 3 x 10 7 K due to the electron ex- 
citation of resonance transitions in the commonest metal 
ions; iii) thermal Bremsstrahlung; and iv) Compton heating- 
cooling. They conclude that thermal instability cannot form 
the quasar line-emitting clouds for the observed pressure 
P ~ 1CP 2 dyne cm~ 2 . On the other hand, magnetic fields 
and intense radiation may provide an important contribu- 
tion to the formation and confinement of the clouds (Rees 
1987; Gongalves et al. 1993a, 1996). The additional heat 
source considered here, AH, together with more conventional 
ones (see Mathews & Ferland 1987; Krolik 1988; Mathews & 
Doane 1990), maintain the broad- line emitting clouds and 
the intercloud medium in pressure equilibrium (Gongalves 
ct al. 1993a, 1996), as required in the two-phase BLR model 
(see Netzer 1990). Within the framework of our model, the 
correlation between continuum and line emission of the BLR 
is also affected by AH as an additional heating/excitation 
source of the BLR, acting together with photoionization by 
the central source. 

The present study extends the stationary equilibrium 
calculations of Gongalves et al. (1993a, 1996, hereafter Pa- 
pers I and II, respectively), in which it was found that AH 
could establish a stable two-phase medium in BLRs, by 
performing time-dependent calculations for the evolution of 
line-emitting clouds from their formation from a 10 7 K in- 
tercloud gas to ~ 10 4 K clouds. We compare the results of 
our models with observational data coming from optical and 
UV AGN monitoring campaigns. 

Section 2 reviews the knowledge of magnetic fields in 
BLRs and the AH mechanisms. The time-dependent models 
for BLR clouds are discussed in Section 3. Model results 
and their comparison with the observations of AGN UV and 
optical lines are given in Section 4. An overall discussion as 
well our main conclusions are stated in Section 5. 



2 ALFVEN HEATING IN THE BLR 
INTERCLOUD MEDIUM 

The classical observational signature of magnetic fields in 
AGN is the radio emission resulting from synchrotron radi- 
ation of relativistic electrons spiralling in a magnetic field. 
Rees (1987), considering the magnetic confinement of BLR 
clouds, has estimated a magnetic field strength of ~ 1 G, 
in the inner parsec of AGN (this value refers to the hot in- 
tercloud gas). Since magnetic field are present in BLRs, we 
also expect the existence of Alfven waves. On disturbing a 
uniform magnetic field, the magnetic tension makes a wave 
propagate along the field lines with the Alfven speed, va 
(for a review see Priest 1994). The perturbations of the mag- 
netic fields could originate in the central region of the AGN, 
which is highly perturbed as we infer from the fast contin- 
uum variability. Therefore, two conditions that are fulfilled 
by the central region of AGN make them favourable sites for 
the generation of Alfven waves: they contain magnetic fields 
and are highly perturbed. 

The diagnosis of the basic physical parameters of BLRs 
is not simple because electron densities are high enough 
for almost all forbidden lines to be collisionally suppressed. 
However, the temperature is estimated via similarity with 



other ionized gases, with temperatures of about 10 4 K. The 
characteristic broad- line widths are ~ 5000 km s _1 , but the 
above temperature implies a line-of-sight velocity dispersion 
of ~ 10 km s _1 ! Therefore, line widths do not reflect thermal 
motions only, but are also due to bulk motions of individ- 
ual broad-line emitting clouds. These clouds are 'usually as- 
sumed, without justification, to have no turbulent velocity' 
(cf. Baldwin 1997). On the other hand, the clouds could not 
be highly turbulent since in this case the effects on the line 
asymmetries would be enormous. The blue and red asym- 
metries of the broad lines are usually explained by infall or 
outflow, but we are still far from a complete explanation of 
these features (see Espey 1997 for a brief review of velocity 
diagnostics in BLRs). Therefore, the line widths reveal the 
existence of bulk motions and, possibly, of turbulence in the 
BLR, thus fulfilling the conditions for the BLR to produce 
Alfven waves. 

Once we admit the existence of Alfven waves in the 
AGN BLRs, the Alfven waves are subject to several damp- 
ing processes. In this work we consider the nonlinear and 
the turbulent damping processes. Assuming spherical geom- 
etry for the clouds, we present below the dependence of the 
Alfven wave heating rate on density and temperature (see 
Papers I and II for a detailed discussion). 

High amplitude waves cause nonlinear mode couplings 
and, in this process, the resulting modes are rapidly damped. 
In general, the waves that result from the damping pro- 
cesses are sound waves; thus, the nonlinear Alfvenic energy 
is converted into thermal energy. An Alfven wave is likely 
to dissipate because of its nonlinear interaction with ei- 
ther a non-uniform ambient field or another Alfven wave 
(Wentzel 1974). The nonlinear interaction of magnetohy- 
drodynamic waves has been treated in detail by Kabu- 
raki & Uchida (1971), Chin & Wentzel (1972) and Uchida 
& Kaburaki (1974). In regions of strong (va > c B ) mag- 
netic field, one Alfven wave can decay into another Alfven 
wave and a sound wave travelling in the opposite direc- 
tion. The resulting Alfven wave has a frequency smaller 
than the original one and it can, in turn, decay into an- 
other lower-frequency Alfven wave plus an acoustic wave. 
The cascade continues until all the Alfvenic energy has been 
converted to acoustic waves that dissipate rapidly. Consid- 
ering this mechanism, in our models, the onset of AH occurs 
for (3 — (B 2 /8tv) /nksT > f3 on > 1, where f3 on is the initial 
value of [3 when AH becomes important (/? > 1 corresponds 
to va/c s > y/6/5 for an ideal 7 = 5/3 gas, where c s is the 
sound speed). Note that our definition of (3 corresponds, in a 
plane parallel geometry, to /3 — ps /p g as , whereas in plasma 
physics, /3 p i = Pgas/pB denotes the "/3 parameter" of the 
plasma. The present convention for /3 has been widely used 
in cooling flow studies (e.g. David & Bregman 1989), and 
has been employed by Friaga et al. (1997) in the study of 
the effects of AH in optical filaments in cooling flows. 

According to Lagage & Cesarsky (1983), the nonlinear 
damping rate is 

where £ = 5 — 10, u) is the characteristic Alfven frequency 
and p(5v 2 ) /(B 2 /8n) the ratio of the energy density of Alfven 
waves to that of the magnetic field. In laboratory plasmas, 
the cut-off frequency of Alfven waves is often comparable to 
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the ion-cyclotron frequency, u) c i = eB/mic, so we can write 
a) as a fraction F of u) C i, i.e., u> = Fu> C i- In addition, there 
are experimental grounds to set F w 0.1 or smaller (Burke, 
Maggs & Morales 1998). 

In terms of AH, e.g. H n i (erg cm' 3 s" 1 ), we have: 
H n i = ($w r n i)/wA, where $ w is the wave flux, T n i is given 
by (1), and va is B/^/Airp. As the perturbation collapses, 
the component of the magnetic field perpendicular to the di- 
rection of compression is amplified. For spherical geometry 
B oc A^ 1 oc p 2 ^ 3 , where A is the cross-sectional area perpen- 
dicular to the magnetic field. As a consequence, va oc p 1 ^ 6 . 
Let the dependence of $ w on p be given by $ w oc A^ 1 oc 
p 2//3 . Taking $ w = p{5v 2 )va, we have p(Sv 2 ) oc p 1 ^ 2 . Since 



oc T ' , for w independent of the density, we obtain: 

1/2 rpl/2 



H nl oc p T 



(2) 



In addition, the nonlinear interaction between outward 
and inward propagating waves results in an energy cas- 
cading process (Tu, Pu & Wei 1984). Following Hollweg 
(1986, 1987), the volumetric heating rate associated with 
the cascade process (e.g. associated with turbulent damping 
of Alfven waves) can be written as 



H tu = p(Sv 2 f /2 /L c . 



(3) 



where L C orr is a measure of the transverse correlation length 
of the magnetic field. Again, for spherical geometry, {5v 2 } oc 

oc B- 1 ' 2 



,1/2-1 

obtain 



p~'~ ' oc p 1//2 , and adopting L, 



oc p 1 ^ 3 , wc 



rr 7/12 
-Htu OC p . 



(4) 



We would like to call the reader's attention to the fact 
that Paper II has a misprint in equations (4) and (6), more 
specifically, in the dependence of $ w on the density (in the 
calculations the correct expression was used). Equations (2) 
and (4) now exhibit the right dependence of $ w on the den- 
sity. As a matter of fact, these heuristically derived density 
dependencies were included in Paper II, since we have dis- 
cused not only AH density dependency due to nonlinear, res- 
onance surface and turbulent Alfvenic dampings, but also a 
range of other density dependencies for "like- Alfvenic heat- 
ings" (see Fig. 9 of Paper II). In the other two papers dis- 
cussing AH in other environments than those of BLRs, the 
correct expressions for H n \ and H tu are found (Friaga et al. 
1997; Gongalves, Jatenco-Pereira & Opher 1998). 



3 TIME-DEPENDENT COOLING 
CONDENSATIONS 

We have investigated the evolution of quasar BLR clouds 
from their formation out of the hot phase of the BLR and 
calculated their ultraviolet and optical signatures within the 
scenario outlined in Section 2. The evolution of the cooling 
clouds is obtained by solving the hydrodynamic equations 
for mass, momentum and energy conservation (see Friaga 
1986, 1993): 



dp 
dt 



dr 



du ^ du 
dt dr 



(r 2 pu) = 

1 dp t _ GM(r) 
p dr r 2 



(5) 



(6) 



dU dU 



(■ 



J_d_ ( Kr 2dT\ 
pr 2 dr V dr J 



(7) 



where u, p, p t ,U,T and k are the gas velocity, density, total 
pressure, the specific internal energy, temperature, and the 
thermal conductivity coefficient. The equation of state , 



y _ 3 k B T wb 
2 fim H p 



(8) 



relates U to T and hb = B 2 /8-k, the internal magnetic en- 
ergy density (&b is the Boltzmann's constant, ran is the 
hydrogen atomic mass and p = 0.62 is the mean molecular 
weight of a fully ionized gas with solar abundances) . 

The total pressure, p t , is the sum of the thermal and 
magnetic pressures, p t = p+p B - The constraints on the mag- 
netic pressure come from considerations on the /3-parameter, 
assuming that the magnetic pressure of the unperturbed 
medium would not exceed its thermal pressure, in accor- 
dance with Mathews & Doane (1990). We assume frozen-in 
fields tangled on scales less than that of the perturbation, 
so that, for spherical collapse, B oc p 2 ^ 3 (pb oc p 4 ^ 3 ). In 
this case the magnetic field can be considered as effectively 
isotropic and we treat the isotropic tangled magnetic field 
as a 7 = 4/3 gas, and so ps — «b/3 = B 2 /24tt. 

The self-gravity of the condensations is taken into ac- 
count through M(r), the mass distribution of the gas. Note, 
however, that the self-gravity is much less important than 
the tidal force due to the BH, as one can see from the com- 
parison between the surface gravity of the perturbation (of 
mass M = 47rp m 7? 3 and radius R = L/2 for the values of p m 
and L defined in the next section) 



g sut = GM/R 2 



8.4 x 10~ 8 cms" 1 



and the tidal acceleration due to a central supermassive 
black hole of mass Mbh = 10 8 Mq at a distance r — 30 
light-days of the perturbation: 

gtid = GM B uR/r 3 = 0.71 cms -1 . 

In the energy equation, the heating term H(p, T) fol- 
lows either equation (2) or equation (4) and it is turned off 
in the case of no heating. In the adopted definition of the 
cooling function A(T), A(T)p 2 , is the cooling rate per unit 
volume. Since there is no ionization equilibrium for temper- 
atures lower than 10 6 K, the ionization state of the gas at 
T < 10 6 K is obtained by solving the time-dependent ion- 
ization equations for all ionic species of H, He, C, N, O, Ne, 
Mg, Si, S, Ar and Fe. The radiative cooling function was 
obtained via the collisional ionization code MEKA (Mewe, 
Gronenschild & van der Oord 1985). The coefficients of col- 
lisional ionization, recombination and charge exchange of 
the ionization equations is also from MEKA. The adopted 
abundances are cosmic abundances usually assumed in AGN 
models (Netzer 1990). If solar abundances are used, they are 
taken from Grevesse & Anders (1989). 

The last term of eq. 7 accounts for the thermal con- 
ductivity. We take the thermal conductivity coefficient, k, 
as a fraction, rj, of the Spitzer (1956) value ks = 4.87 x 
10~ 7 erg enf 1 K _1 , with t) = 10~ 4 . The thermal conduc- 
tivity should be inhibited, otherwise it would prevent the 
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development of the thermal instability. The thermal conduc- 
tivity would erase any steep temperature gradient arising in 
the medium, thus heating the cooler gas of BLRs to the 
higher temperatures of the surrounding hot gas. Therefore, 
the presence of colder clouds embedded in a hot gas requires 
some reduction of the thermal conduction. A similar inhibi- 
tion of thermal conduction seems to be needed also in cool- 
ing flows in clusters of galaxies: the classic (Spitzer) thermal 
conductivity would drive a heat flow towards the cooling 
flow, and, as a consequence, there should be no cooling flow! 
The magnetic fields could be at the origin of the inhibition of 
thermal conductivity. In the presence of magnetic fields, the 
heat conduction is suppressed in the direction perpendicular 
to the magnetic field, taking place substantially only along 
magnetic field lines. If the magnetic field is very tangled, 
the heat conduction can be reduced by a large factor. Other 
processes of inhibition of thermal conductivity besides mag- 
netic field entanglement are possible. For instance, in cooling 
flows in galactic clusters, electromagnetic instabilities driven 
by temperature gradients (or electric currents in other situ- 
ations) can also cause this inhibition, even for non-tangled 
field lines (Pistinner, Levinson & Eichler 1996). 

The spherically symmetric hydrodynamical equations 
are solved using a finite-difference, implicit code based on 
Cloutman (1980). The code is run in the Lagrangian mode 
and the grid points are initially spaced logarithmically, with 
a grid of 100 cells. The artificial viscosity for the treatment of 
the shocks follows the formulation of Tsharnutter & Winkler 
(1979) based on the Navier-Stokes equation. In contrast with 
the von Neumann-Richtmyer artificial viscosity (Richtmyer 
& Morton 1967), this form of artificial viscosity vanishes 
for homologous contraction. Thermal conduction is included 
by operator-splitting using a Crank-Nicholson scheme. The 
outer boundary conditions on pressure and density are de- 
rived by including an outer fictitious cell, the density and 
pressure in which are obtained from extrapolation of power 
laws over the radius fitted to the outermost real cells. 

The initial density perturbations are characterized by 
an amplitude, A, and a length scale, L, following 



Sp 
P 



Asin(x) 



2-wr 



(9) 



We have also assumed that the perturbations are isobaric 
and nonlinear (^4 =1), avoiding the uncertainties concerning 
the linear development of thermally unstable modes in the 
presence of a gravitational field (David, Bregman & Seab 
1988). The unperturbed hot-phase intercloud medium is 
characterized by a total hydrogen density n m = 5 x 10 6 cm~ 3 
and a temperature T m = 10 7 K, appropriate for the in- 
ner regions of an AGN, ~ 10 18 cm from the central source, 
where lies the BLR. If the gas cools isobarically, when it 
reaches T = 10 4 K, the resulting central electronic den- 
sity would be 6 x 10 9 cm~ 3 , consistent with the ab- 
sence of forbidden line emission in AGN BLRs. We fixed 
0o = B 2 /{8 TrnksTm) = 0.03 in accordance with Mathews 
& Doane (1990) who pointed out that this quantity should 
be between 0.003 and 0.03, in a similar calculation with no 
Alfvenic heating. The length scale of the perturbation was 
fixed at 5 x 10 16 cm. The corresponding column density is 
Ah = 3.2 x 10 23 cm" 2 . 

Here we investigate several representative models with 
magnetic fields, one with no AH and the others with AH 



due to two damping mechanisms (nonlinear and turbulent) 
and at different efficiencies (see Table 1). Model I allows 
for the presence of a magnetic field only through a magnetic 
pressure term in the equation of motion. In models II and III, 
the energy equation includes AH as an additional heating 
term. We considered nonlinear and turbulent Alfven wave 
heating, from the above equations (eqs. (2) and (4)), and 
these heating terms have the forms 



tf nl = C#o(n/n )- 1/2 (T/To) 1/2 , 
for nonlinear heating, and 
H tu = CH (n/n ) 7/12 , 



(10) 



(11) 



for turbulent heating, where £ is the efficiency of Alfven 
heating. The choice of the normalization is the following: 1) 
To = 10 5 K (the temperature at which optical emission be- 
gins to be important); 2) no = 13.59 x (5 x 10 6 ) cm -3 , where 
13.59 would be the compression factor from the unperturbed 
n m = 5 x 10 6 cm -3 and T m = 10 7 K initial state to T = T , 
under isobaric conditions for a plane parallel geometry; and 
3) Ho = (3/2) noksT /t co \ erg cm~ 3 s _1 , where t co \ is the 
collapse time (defined below) for the cloud with no Alfvenic 
heating. 



4 MODELLING RESULTS AND DATA 
COMPARISON 

We investigate models with magnetic field and AH (excepted 
model I in which there is no AH) using H n i and H tu with 
the efficiency of the conversion of Alfven wave energy in 
Alfvenic heating, £, equivalent to 10% and 30%. In model I 
the magnetic effect is considered only as magnetic pressure 
in the equation of motion. All the other models contain AH 
as an additional heating term in the energy equation. The 
characteristics of our models are summarized in Table 1: 
model identification (column 1); the Alfvenic heating con- 
sidered in each model (column 2) ; the efficiency £ of the AH 
(column 3); and the total collapse time tt.coi, i-e. the time 
since the beginning of the calculations until the time when 
the innermost cell reaches 8 x 10 3 K (column 4). 

First of all, note that AH can maintain stable equilib- 
rium between the optical line-emitting cloud and the inter- 
cloud medium. In fact, Friaga et al. (1997), studying the 
Alfven heating/excitation of the optical filaments in cooling 
flows, have found that one of their models (with nonlinear 
heating and ( — 100%) reached thermal stability, showing 
that AH can suppress strong growth of thermal instabili- 
ties in cooling flows. This suggests that AH could induce a 
stable two-phase equilibrium in other astrophysical environ- 
ments, such as the BLRs of AGN, which were investigated 
with simple time-independent calculations in Papers I and 
II. Therefore, we take the initial conditions of our simula- 
tions from Paper II. As we see from Table 1, an increase in 
the AH efficiency from 10% to 30%, leads to longer total 
collapse times: from 10.10 to 10.60 xlO 7 s for the models 
with H n i and from 9.89 to 9.93 xlO 7 s for those with H tu - 
The same trend of collapse time with £ was found by Friaga 
et al. (1997), thus confirming that AH inhibits development 
of thermal instabilities. 

Figure 1 -3 show the temperature, thermal pressure and 
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Table 1. Properties of the Models 



Model AH C t t , co i 

10 7 s 

I none - 9.88 

Ilnl nonlinear 10% 10.10 

Iltu turbulent 10% 9.89 

Illnl nonlinear 30% 10.60 

Illtu turbulent 30% 9.93 




10 1B 10 lfl 

model Illtu 



10' r 




10 15 10 ,s 

Distance (cm) 



Figure 1. Profiles of (a) temperature for models I(no AH), 
IInl(H nl , C = 0.1), IItu(H tu , C = 0.1), IIInl(H nl , C = 0.3) and 
IIItu(_f/tu, f = 0.3) for several evolution times. 



P (magnetic to thermal pressure ratio) profiles of our mod- 
els. These profiles are shown from the beginning of the cal- 
culations, throughout the optical phase (T < 10 s K), until 
9.8 x 10 7 s (very close to the total collapse time of model 
I, 9.88 x 10 7 s). The influence of the AH on the collapse of 
the perturbed gas can be well understood following the pro- 
files of Figure 1. For model I (with no AH), nearly all the 
gas contained in the perturbation reaches ~ 10 4 K in the 
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total collapse time of 9.88 x 10 7 s. All the models with AH 
have total collapse times longer than this. Comparing the 
results for nonlinear and turbulent AH models, we can see 
that nonlinear AH is more effective in inhibiting the devel- 
opment of the collapse than turbulent AH. Note that model 
Illnl (nonlinear AH, with £ = 30%) has the highest total 
collapse time, and therefore the smallest amount of "cool" 
gas when the evolution time reaches 9.8 x 10 7 s. 

As we see from the thermal pressure profiles of Fig- 
ure 2, as compared to Figure 1, the drop in the temperature 
causes a reduction in the thermal pressure. This behaviour 
is a consequence of the fact that the sound crossing time 
is much longer than the cooling time, so there was no time 
for repressurizing. The (3 profiles (Figure 3) show that, at 
the beginning of the calculations, the magnetic pressure is 
much smaller than the thermal one; but, as the perturbation 
evolves, the magnetic pressure becomes more important, as 
a result of the drop in the temperature of the perturbed gas, 
P becomes > 1, and any (slight) compression of the gas will 
contribute to the increase in f3 (assuming frozen-in fields, 
the magnetic field increases according to B oc p 2 ^ for a 
spherical geometry). 

When the temperature falls to ~ 5 x 10 K, the gas 
begins to emit optical and UV lines. Before this optical line 
emitting stage, the gas cooling occurs via EUV and X-rays 
lines. We have calculated the emissivity of low- ionization 
optical lines as well as that of high-ionization UV lines. 
The results of our models may be compared to the observa- 
tional data from monitoring campaigns, such as those of the 
Seyfert 1 galaxies Fairall 9 (with IUE, Rodrfguez-Pascual et 
al. 1997; in the optical range, Santos-Lleo et al. 1997) and 
NGC 5548 (Dumont et al. 1998). 

We set the beginning of the phase of optical-UV line 
emission when the temperature of the innermost cell (that 
is, the first cell to cool) falls below 5 x 10 5 K. Figure 4 
shows the evolution during this phase of the luminosity of 
the low-ionization line H/3, and of the high-ionization lines 
Hen 1640 A, Civ 1549 A, Nv 1240 A and Ovi 1034 A, for the 
models I, Ilnl and Illnl. The variability pattern of the three 
models is very similar for all the lines studied. In addition, 
AH enhances significantly, albeit not strongly, the broad line 
emission. We should note that the amplitude of variation of 
the H/3 line is smaller than that of the high-ionization ions, 
and that the time-scale of its variability is much longer. 

With respect to the variability of the collisionally ex- 
cited high-ionization UV lines, their luminosities show well- 
defined maxima occurring ~ 10-25 days after the appear- 
ance of the first gas with T < 5 x 10 5 K. The variability 
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Figure 2. Profiles of (b) thermal pressure for models I(no AH), 
IInl(H n i, C = 0.1), IItu(if t u, C = 0.1), IIInl(J? n i, C = 0.3) and 
IIItu(/ftu, ( = 0.3) for several evolution times. 



Figure 3. Profiles of (c) magnetic to thermal pressures ratio 
for models I(no AH), IInl(if n i, £ = 0.1), IItu(H tu , C = 0.1), 
IIInl(ff n i, ( = 0.3) and IIItu(// tu , ( = 0.3) for several evolution 
times. 



of the lines Civ 1549 A, Nv 1240 A and Ovi 1034 A is 
faster and of larger amplitude than that of the Hell 1640 A 
line. For model Illnl, the time interval Ati/2 for the line 
luminosity to decrease from its maximum value to half- 
maximum value is 46.0 d, 24.2 d, and 31.9 d, for the lines 
Civ 1549 A, Nv 1240 A and Ovi 1034 A, respectively, 
whereas for Hell 1640 A, At 1 / 2 = 152 d. Note that, in ad- 
dition to increase the luminosities of the lines (during the 
period from the time of peak luminosity and t = 100 d, the 
increase is 30-55% for Hell 1640 A, 60-70% for Civ 1549 A, 
48 - 52% for Nv 1240 A, and 51 - 60% for Ovi 1034 A) the 
AH makes the luminosity peaks more prominent. 

The evolution of the H/3 luminosity contrasts with the 
evolution of the high-ionization UV lines in that it is almost 
constant within an elapsed time of about 300 days. It shows 
a broad maximum at 150-200 days and varies within only 
10% between 70 and 320 days. The fact that the amplitude 
of variation of the H/3 line is smaller than that of the high- 
ionization ions, and that the time-scale of its variability is 
much longer is in a first order agreement with observations. 
For Fairall 9, the H/3 line shows a smoother pattern of vari- 
ability than Hen 1640 A, with a time-lag r d = 23.0 d and a 



ratio of maximum to minimum fluxes i? max = 1.55, whereas 
for Hen 1640 A, r d = 4.2 d and i? max = 3.25. However, a 
closer comparison of the luminosity and time-lags between 
Balmer and high-ionization lines reveals that the predicted 
luminosity of H/3 is too low and its variability too slow. For 
model Illnl, at t — 150 d, even close to the maximum of the 
H/3 line, and far from the maximum of the Civ 1549 A line, 
the ratio Civ 1549 A/H/3 is 24.7, whereas the typical value 
for AGN BLRs is ~ 5. AH does not seem to account for 
the Balmer recombination lines, which could be explained 
by photoionization, but it makes a significant contribution, 
both to the energetics and to the variability of the high- 
ionization lines. 



5 DISCUSSION AND CONCLUSIONS 

In our model, the values of the time of the luminosity peak 
= 29.7 d for Civ 1549 A, £ max = 21.9 d for 
Nv 1240 A, and £ max = 21.9 d for Hell 1640 A, for model Il- 
lnl — reproduce the observed sequence of time-lags in AGN 
BLRs, with the time lags for Nv 1240 A and Hell 1640 A 
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Figure 4. Luminosity variability of the line H/3, as well as those of 
the UV high-ionization lines Hen 1640 A, Civ 1549 A, Nv 1240 A 
and Ovi 1034 A, for models I (dotted lines), Ilnl (dashed lines) 
and Illnl (solid lines). The time is counted from the beginning of 
the optical phase (i.e. when T becomes < 5 X 10 5 K). 



being of comparable magnitude and shorter than those of 
Civ 1549 A ( our models predicts an even shorter variability 
time scale for Ovi 1034 A: i max = 11.7 d for model Illnl). 
The usual interpretation is that the Hell-Nv emission region 
is located at a shorter distance from the central continuum 
source than the Civ 1549 A region — in the case of NGC 
5548, at 6 and 10 light-days, respectively — (Dumont et al. 
1998). The present results suggest that at least part of this 
trend could be due to different evolutionary time scales for 
each line within one single cloud or cloud ensemble at one 
given position with respect to the central source, instead of 
different clouds giving rise predominantly to different lines 
at differing locations. 

We can interpret our results in the light of the discus- 
sions by Dumont et al. (1998), which, although give support 
to scenario in which the BLR consists at least of two different 
regions, in which low- and high-ionization lines are predom- 
inantly produced, also point out some inconsistencies of the 
photoionization models based in this scenario with respect 
to the energy budget, the line ratios and the line variability. 
The variability represents the greater challenge for the strat- 
ified model for the formation of the high- and low-ionization 
lines, in which the more ionized lines are formed closer to 
the center (CIV 1549A and He II 1640A) and the less ion- 



ized further out (Collin-Souffrin et al. 1986; Collin-Souffrin, 
Hameury & Joly 1988). 

It is possible that the standard stratified model suffers 
of an excessive number of emission regions. For NGC 5548, 
Dumont et al. (1998) distinguishes four zones: the Mgll- 
CIII] emission zone, at ~ 30 — 50 lt-days; the Balmer lines- 
Hel emission zone, at ~ 12 — 20 lt-days; the Lycv-CIV emis- 
sion zone, at 8 — 12 lt-days; the NV-Hell emission zone, at 
4 lt-days. Our results also allow for several locations for the 
formation of the broad lines, with the high-ionization lines 
being produced closer to the central continuum source, and 
having an important contribution from AH, and the Balmer 
lines being produced at a more distant location, and pow- 
ered mostly by photoionization, which would be responsi- 
ble for most of its variability. Of course, some fraction of 
the broad-line spectrum could be non-variable, and the AH- 
powered clouds could contribute to a slowly-variable part of 
the Balmer lines emission. The AH mechanism would be one 
more ingredient in the constitution of the time-lags. The de- 
creasing values of t max for the lines Civ 1549 A, Nv 1240 A 
and Ovi 1034 A shows that different time lags do not nec- 
essarily indicate different line formation regions. A strati- 
fied model simpler than the standard one could consider the 
formation of the low-ionization lines in clouds far from the 
center accommodate the production of the highest ioniza- 
tion lines (from Civ 1549 A to Ovi 1034 A) in the same 
location very close to the center (e.g. in the case of NGC 
5548, in the few inner lt-day), with part of delay of the lines 
with respect to fluctuations of the continuum arising from 
evolutionary effects within the clouds, and heated/excited 
by non-radiative processes, in our case, by AH. 

The use of Alfven heating in our models is not aimed 
at solving all the problems pointed out by Dumont et al 
(1998): the 'energy budget problem', the 'line ratio prob- 
lem' and the 'line variation problem'. For example, AH is 
no solution to the 'energy budget problem', because of the 
small contribution of AH to Ha emission (and to the Balmer 
series). However, AH alleviates the 'line variation problem' 
because it reproduces the trend of time lags with different 
ionization lines, with no need for an excessive stratification 
of clouds emitting different lines. Although AH does not ac- 
count for the Balmer recombination lines, which could be 
explained by photoionization, it makes a significant contri- 
bution both to the energetics and to the variability of the 
high-ionization lines. What we explored here is the possi- 
bility that several species of high-ionization lines, are pro- 
duced at the same location, however, being not only pho- 
toionized but also heated/excited by, in the present case, 
Alfvenic heating. 

The lag between continuum and line variability could 
arise from some process associated with maxima in the con- 
tinuum emission initiating the collapse. A possible mecha- 
nism for triggering the collapse of the perturbation is Comp- 
ton cooling. Compton cooling resulting from luminous QSOs 
in clusters of galaxies has been shown to be able to induce 
or sustaining cooling flows (Fabian & Crawford 1990; Atuko 
et al. 1998). Compton cooling can also be important in ac- 
cretion disks (Esin 1997). Mathews and Ferland (1987) find 
that the Compton temperatures appropriate to quasar spec- 
tra are in the range 2.8 x 10 6 < T c < 4.6 x 10 7 K. If the 
source radiation is soft (Tc < T m , the temperature of the 
intercloud medium), burst in the central continuum could 
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induce cooling on time scales close to the Compton cooling 
time: 



tc = 



lnk B T m 



fc B CT T _^bol 
" m c c' 2 i-nFi 2 



n(T m - T c ) 



3 m e c 2 R 2 Tc.-t 

I (XT ivbol -*m 



where 7? is the distance to the source and Lboi is its bolo- 
metric luminosity (and the other symbols have their usual 
meaning). For typical values of R and Lboi in AGN BLRs, 



t c = 4.06 



(R/30 light days) 2 



(L bol /10«ergs-i) (1 tJ 



days, 



so it is feasible to induce collapse of perturbations on the 
required short time scales. 

As mentioned in section 2, high amplitude waves cause 
nonlinear couplings and the nonlinear Alfvenic energy is, in 
general, converted into thermal energy. The nonlinear damp- 
ing rate (eq. (1)) is a function of Q, which can be written in 
terms of the ion-cyclotron frequency, uj — Fuj c i. Wc assume 
F = 0.1, which is consistent with laboratory experiments on 
heat transport in magnetized plasmas. In a electron beam 
experiment (Burke et al. 1998), classical conduction is ob- 
served for the first 2 ms since beam turn-on, after which 
fluctuations with frequencies ~ O.lWci spontaneously grow, 
and at 3 ms, large amplitude fluctuations with frequencies 
0.02cj c i appear. Following equation (1), the damping length 
Idamp = va/Y can be written as 



J damp 



)- 1 F- 1 n H 1/2 /3 1 / 2 cm. 



For the onset of the AH (/? = 1, n H = 5 x 10 6 cm 3 ), as- 
suming i = 5, P (Sv 2 )/(B 2 /8ty) = 1, and F = 0.1, L damp = 
8.7 x 10 4 cm, much smaller compared with the initial dimen- 
sion of the system (5 x 10 16 cm) and the radius where the 
Alfvenic heating is acting (4.2 x 10 15 cm for the model Illnl). 
This implies that the waves are damped very rapidly and so 
must be generated locally in the cloud. We have evidence for 
generation in situ of Alfven waves from observations of out- 
ward and inward fluctuations propagation in the solar wind 
(Coleman 1968; Belcher & Davis 1971). Roberts et al. (1987) 
using data from Helios and Voyager spacecraft, investigated 
the origin and evolution of low-frequency planetary fluctu- 
ations from 0.3 to 20AU. They concluded that the outward 
traveling fluctuations are predominantly generated by the 
Sun, but that in situ turbulence, most likely due to stream 
shear, generates fluctuations with both inward and outward 
senses of correlation. This mixed state provides stronger ev- 
idence for in situ generation. We assume that the above sce- 
nario can be applied to the BLR clouds allowing the in situ 
generation of Alfven waves. 

It would be interesting to compare the results of the 
present model with other models assuming local heating 
sources. Besides photoionization model by central contin- 
uum, local emission models have also been proposed in or- 
der to explain the Ha luminosities and the major line ratios 
of cluster nebulae (Voit et al. 1994). The major difference 
between the work of Voit et al. (1994) and the present one 
is that we are dealing with an evolutionary model where the 
heating source is non-radiative (AH heating) and the total 
pressure includes both thermal and magnetic pressures. On 
the other hand, in the model of Voit et al. (1994) it is as- 
sumed a pressure and heating-cooling equilibrium, while our 



model is time-dependent and does not assume any equilib- 
rium. Moreover, Voit et al. (1994) present the results of a 
photoionization code (CLOUDY), considering the soft X- 
ray/EUV radiation field of a cooling flow in which cold 
clouds are embedded, while we restricted ourselves to purely 
collisional emissivity calculations. 

Our model is not aimed at reproducing the continuum 
emission of any particular AGN. Instead is a tool to in- 
vestigate separately from photonization effects the role of 
magnetic fields (magnetic pressure and Alfven heating) as a 
mechanism contributing to the delay between the continuum 
variability and the line variability. A complete model for the 
BLRs would include photoionization, shocks and magnetic 
effects. The focus of the present work is on the third ingre- 
dient. 

We should also note that even non-standard radiative 
heating mechanisms deserve proper consideration. If the hot 
intercloud medium is turbulent, mixing layers are expected 
to develop around the colder BLR clouds. In the conditions 
of cooling flows, Begelman & Fabian (1990) have shown that 
EUV radiation from mixing layers could power the line emis- 
sion of optical filaments in cooling flows. It is possible that, 
also in BLRs, the radiation from mixing layers around the 
BLR clouds could constitute an additional, local, photoion- 
ization source, besides the central continuum. 

Other issue is the stability of the BLR cold clouds. 
Several external forces are expected to be acting on (and 
against) the BLR clouds in the harsh environment con- 
stituted by the intercloud medium — high temperatures 
(T ~ 10 7 ~ 8 K), high relative velocities (up to « 10 4 km 
s _1 ), and, probably, turbulence. Significant motions of the 
BLR clouds relative to the intercloud medium could lead to 
disruption of the clouds, and radiation forces could ablate 
the protoclouds before they contribute fully to the broad 
line emission (Mathews & Doane 1990). Also with respect 
to this topic, magnetic fields have a role to play, providing 
supporting pressure and confinement of the clouds (Rees 
1987; Gongalves et al. 1993, 1996). In addition, in the mo- 
tion of the BLR clouds through the whole of BLR (with the 
accompanying effects on the line profiles) magnetic fields 
could have a dynamical role, as it seems to the case in cool- 
ing flows (Gongalves & Friaga 1999). Finally, in order that 
the density reaches values high enough to account for the ab- 
sence of forbidden lines in AGN BLRs, some process should 
remove the magnetic field in the late evolution of the per- 
turbations, when pressure equilibrium has been achieved, 
otherwise the magnetic field would prevent further compres- 
sion of the BLR cloud. One mechanism that was shown to 
be effective in optical filaments in cooling flows is magnetic 
reconnection (MR) (Jafelice & Friaga 1996, Friaga & Jafelice 
1999). Therefore, it would be interesting to investigate MR 
in BLR clouds, both as a heating source and as mechanism 
for removing magnetic support. 
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